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Abstract 

We investigate low-energy cosmic-ray antiprotons (p's) arising from the frag- 
mentation of quarks and gluons emitted from evaporating primordial black 
holes (PBHs). To calculate the local interstellar flux of these p's, their prop- 
agation in the Galaxy is described by a 3-D Monte Carlo simulation based on 
the diffusion model. This flux is used with recent observations to derive new 
upper limits on (i) the local PBH explosion rate 7^ < 1.7 x 10^^ pc^^yr"^, 
(ii) the fraction of the Universe's mass going into PBHs with particular mass, 

and (iii) the average density of PBHs in the Universe. 
PACS numbers: 97.60.Lf, 96.40.-z, 95.85.Ry, 98.35.Gi 
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From a standpoint of cosmological understanding, it would be of great value to confirm 
the existence or nonexistence of primordial black holes (PBHs), which may have formed 
in the early Universe via initial density fluctuations, phase transitions, or the collapse of 
cosmic strings (for a review, see Ref. Hawking first showed that black holes (BHs) 
emit particles and evaporate by quantum effects, noting that PBHs are the only ones with 
a mass small enough for the quantum emission rate to be significant, possibly yielding an 
observable effect. For example, the hard 7-rays from small enough PBHs may contribute 
to the diffuse 7-ray background spectrum, though no distinct signature has been observed. 
Thus, this leads to an upper limit (U.L.) on the average density of PBHs in the Universe, 
i.e., the ratio of their density to the critical density of the Universe, I^pbh, must be ^ 10~^ 

i- 

Despite such a stringent limit, their signature could still appear in the spectrum of 
cosmic-ray antiprotons (p's) [Q, which are generally considered to be secondary products 
from interactions of cosmic-ray protons (p's) with the interstellar medium. This possibility 
arises because, although the kinematics of secondary p production by such p's should lead 
to a steep drop in the resultant p flux at kinetic energies less than 2 GeV, the expected 
flux of p's from PBHs (PBH-p's) has contrastingly been shown to increase with decreasing 
kinetic energy down to ~ 0.2 GeV [Q; thus providing a distinct signature below 1 GeV. 
Hence, searches for such low-energy cosmic-ray p's could lead to a novel constraint on the 
density of PBHs, or more importantly, demonstrate their existence. With this in mind, and 
spurred by the recent detection of cosmic-ray p's with kinetic energies less than 0.5 GeV, 
being accomplished during a 13-h balloon flight (BESS '93) 0, we present a new method 
describing the propagation of p's in the Galaxy; thereby obtaining the most accurate-to-date 
spectrum of local interstellar PBH-p flux. We also use these results with observed data 
to derive a new U.L. on the density of PBHs which are expiring near the Solar system, 
after which we discuss its cosmological aspects assuming their formation via initial density 
fluctuations. 

The source spectrum of PBH-p's is determined using general properties of BH evaporation 
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Briefly, an uncharged, non- rotating BH with mass M emits particles with spin s and 
total energy between (Q, Q + dQ) at a rate 



dN _ TsdQ 
~dt ~ 27rh 



exp 



Q 



\2s 



-1 



(1) 



per degree of particle freedom, where T is the BH temperature (= hc^ /SnGMk = 1.06 x 
10^^(M/g)^^ GeV, with k = 1), and Tg is the dimensionless absorption probability for the 
emitted species. Considering all species, the corresponding mass loss rate can be expressed 
as 

^ = -5.34 X lO^V(M) (^^) ' g s-\ (2) 

where f{M), a function of the number of emitted species, is normalized to unity for large 
M (> 10^^ g) and increases with decreasing M. From Eq. (|^), it follows that M^.. ~ 
5.3 X 10^''(t„/16 Gyr)^/^ g, or 3 x 10~-'^^(t„/16 Gyr)-*^/^ in units of solar mass, where M* is 
the initial mass of a PBH expiring today; i.e., its initial lifetime equals the present age of 
the Universe tu, being taken here as 16 Gyr |^. PBHs with initial mass Mj < M^, should 
have completely evaporated by now, while those with Mj slightly larger than have an 
extremely high present temperature and will soon expire by "explosion." 

Under the assumptions that (i) Eq. (|l]) holds for each emitted species of quarks and 
gluons and (ii) all PBH-p's are their fragments, which is consistent with observed e^e~ 
annihilations, we calculate the PBH-p source spectrum per unit volume via the following 
three-step procedure: (1) The JETSET 7.4 Monte Carlo simulation code is used to obtain 
the fragmentation function dgp\Ep,Q)/dEp describing the fragmentation of each emitted 
species j with total energy Q into p's with total energy Ep; (2) The p emission spectrum 
from a BH at present temperature T, d^p{Ep,T)/dEp, is calculated by convolving Eq. (|I]) 
with dgp\Ep, Q)/dEp for Q > Ep and summing over all j and their degrees of freedom; and 
(3) The expected PBH-p source spectrum per unit volume is obtained by convolving the 
PBH present temperature distribution dn/dT with d^p{Ep,T)/dEp for T ^ 0.1 GeV, where 
n is the number of PBHs per unit volume. Under Eq. (|^), dn/dT versus T 0.1 GeV) is 



roughly a power-law function (oc T~^'^), with its normalization being solely determined by 
the value of the initial mass spectrum dn/dMi at Mj = M^.. 

We checked the reliability of our use of JETSET 7.4 |^ via data from e^e^ colliders, 
i.e., observed p spectra obtained from the fragmentation of quarks and gluons at various Q 
relevant to evaporating PBHs. First, we verified that the predominant contribution to PBH- 
p's occurs aX Q = 1-5 GeV, after which the simulated p spectra from quark fragmentation 



were compared with DASP data at Q = 1.8-2.5 GeV |ig and ARGUS data at Q = 4.99 GeV 
10| , while those from gluon fragmentation with ARGUS data from the direct decay of T(IS'), 
which mainly proceeds through three gluons with an average energy of mT{i5)/3 = 3.2 GeV 
10| . Resultant comparisons showed good agreement (data not shown). 

Since the propagation of PBH-p's in the Galaxy has thus far only been roughly treated 
we calculate the local interstellar PBH-p flux from the PBH-p source spectrum utilizing 
a 3-D Monte Carlo simulation code based on the diffusion model; an expansion of the 1- 
D simulation by Owens and Jokipii [|1T|. We assume that the PBH spatial distribution is 
proportional to the mass density distribution within the Galactic halo, i.e., oc (1 + (r/r^)^)^^ 
|T^, where r is a distance from the Galactic center and is the core radius of 7.8 kpc. As 



shown later, the local PBH-p flux arises only from nearly expired PBHs existing within a 
few kpc away from the Solar system, whose location is at r = = 8.5 kpc. Accordingly, the 
local PBH-p flux can be calculated via the three-step procedure by simply using dn{Mi = 
M^,,r = rQ)/dMi, which is parametrized by introducing an unknown parameter e^,: 

'^''iM, = M,,r = rr,)=s.^, (3) 



dMi' ' Ml 

where phQ is the local density of halo dark matter (~ 0.3 GeV cm"^ [111])- Under Eq. (^, 
represents the ratio of the density of PBHs with Mj = to phQ- 

Details regarding the diffusive propagation of p's will be described elsewhere [|14|. Briefly, 



however, each p ejected from the fragmentation of quarks and gluons is assigned to an initial 



position xq and energy Eq, after which it travels x to x + uVS-DAt + {'VD)At in each 
subsequent time step At, where u\/6DAt represents the effect of isotropic diffusion during 



At, u is a unit vector with random direction, D = -D(x) is the diffusion coefficient, and 
{'VD)At expresses the anisotropy of diffusion caused by the spatial gradient of D. We 
do not consider the effect of convection due to galactic wind, because Webber et al. |T5[ 
showed that it does not significantly affect cosmic-ray propagation. Regarding subsequent 
collisions and energy loss, the interstellar medium (ISM) is considered to consist of 90% 
hydrogen and 10% helium atoms [l^, and its number density distribution is modeled as 



natom = 1.1 exp(— 2;/100 pc) atoms/cm^, where z is the perpendicular distance from the 
galactic plane {z = 0). While propagating through the ISM, p's lose energy by ionization 
and are lost by annihilation. Those p's passing near the Solar system 25 pc) are included 
in the ffux calculation. 



The Galaxy is modeled as a cylindrical diffusing halo with a diameter of 40 kpc and 
halo thickness of 2h. Free escape is assumed to occur at the boundaries. As normally done, 
we parametrize the diffusion coefficient as D = Dq{z)Di{R), where R is the rigidity of the 
particle. We apply three halo models for Dq{z) (Fig. |l](a) shows Do{z)~^): two with constant 
Dq but different values of h (models I and II), and one in which Dq is dependent on z (model 
III). Assuming Dq is constant, Dq/H and Di{R) are determined |T3[ as (8±1.6) x 10^ cm s~^ 
and ~ (i?/GV)°'^, respectively, which fits the secondary to primary ratios of cosmic-ray 
nuclei, e.g., the ratios of boron to carbon (B/C) and sub- iron to iron (sub-Fe/Fe). The value 
of h is also constrained from 2 to 4 kpc in order to fit the radioactive secondary to stable 



secondary ratios, e.g., the ratios of radioactive beryllium to stable beryllium (^°Be/^Be). 
Thus, h = 2 and 4 kpc in models I and II, respectively. Model III also fits the above ratios 
at the same Di{R). 

The simulated distribution of PBHs contributing to the integrated p flux F near the 
Solar system, obtained by assuming e^, = 1.0 x 10~^, is shown for halo models I-III in 
Fig. |l|(b) and (c) in the z- and r-direction, respectively. Note that all models show the same 
contribution at z = (galactic plane) as they are constrained to reproduce the observed 
secondary to primary ratios which originate from secondary production occurring there. In 
addition, only PBHs within a few kpc away from the Solar system contribute a substantial 



flux. When considering these results along with the fact that the mass density within this 
region is relatively constrained |13[ , this indicates that the Solar system p flux is only slightly 



dependent on the assumed mass density distribution within the Galactic halo [|1^]. The p 
mean confinement time is also calculated by the code as 4.0, 7.9, and 8.0 xlO^ yr for models 
I-llI, resulting in F = 0.78, 1.53, and 1.55 x 10^^ m~^s~^sr~^, respectively. 

As we have shown that the local PBH-p flux can only be due to contributions from PBHs 
that are close to explosion and exist within a few kpc away from the Solar system, data from 
searches in which low-energy cosmic-ray p's are detected can directly constrain (or possibly 
reveal) the PBH explosion rate averaged over this local region, termed here as the local PBH 
explosion rate 71, i.e., 

dn 

n = -j-{Ti = tu,r = r©) 
dTi 

dn , , , , 

(M, = M„r = ro)x-^, (4) 



dMC ' ^' 3t„ 

where (= tu{Mi/ M^Y at Mj M*) is initial lifetime of PBHs with initial mass Mj. 
Figure ^(a) shows simulated local interstellar PBH-p flux using model III with two possible 
values of TZ, i.e., 1 x 10~^ and 1 x 10^'^ pc^'^yr^^. Then, to compare this p flux with 
observational data, we converted it into the p/p ratio at the top of the atmosphere (TOA) by 
assuming that (i) the interstellar p flux Jp = 1.5x 10^/3-^(i?/GV)-2-74 in-2g-ig^-iQgV-^ [JT^, 
where (3 is the ratio of particle velocity to the velocity of light, and (ii) the corresponding flux 
of both p's and p's at TOA can be similarly calculated by accounting for solar modulation 
using the force-fleld approximation [|l^ with modulation parameter = 650 MV, which 
corresponds to that obtained from the BESS '93 flight [Q. Figure 0(b) compares the resultant 
p/p ratio for PBH-p's with the expected p/p ratio for secondary p's and observational 



data p,|T9|-Q. Note that no distinct signature of evaporating PBHs is apparent, i.e., at 
kinetic energies below 1 GeV, the data shows no tendency to reach a constant p/p ratio. As 
no signature exists, statistical analysis of recent observations leads to the following U.L. 
on IZ with 90% confldence level (C.L.): 

7^ < 1.7 X 10^2 pc"V~\ (5) 



which is almost eight orders of magnitude more stringent than the present U.L. on the rate 
of 50-TeV 7-ray bursts {TZ < 8.5 x 10^ pc^'^yr"^ P^)? ^"^^ the practical sensitivity for 100- 



MeV 7-ray bursts {TZ ~ 10^ pc '^yr ^ P^)- Note that the observed diffuse 7-ray background 



spectrum, which contains contributions from all evaporating PBHs throughout the Universe, 
does not directly constrain TZ. 

Equations (H)-® give the following U.L. on with 90% C.L.: 



e, < 2.8 X 10"^ 



PhG \ V ^ 



\ 4/3 



.0.3GeVcm-3y Vl^Gyrj ' 

which can be used to derive an U.L. on the fraction of the Universe's mass going into 
PBHs with mass M^,, i.e., /3(M^,) If PBHs are assumed to have formed via initial 
density fluctuations, /3(M=k) should be closely related to their amplitude on a scale of M*. 
Thus, under this assumption, and for a fiat Friedmann universe, Eq. (^) gives (3{M^) ~ 
1 X 10^^^£:*(i7/i/0.1), where Qh is the ratio of the density of halo dark matter to the critical 
density of the Universe. Finally, assuming Qh = 0.1 f^, Eq. (H) leads to an U.L. of P{M^) < 
3 X 10-26. 

Further, by assuming that such initial density fluctuations are scale-invariant, the initial 
mass spectrum of PBHs should have a power-law form, i.e., dn/dMi cx M"" |l|, where 
a = 5/2 if the Universe was radiation-dominated when PBHs formed. Using this initial mass 
spectrum normalized by Eq. at Mi = M*, we integrate Mj x dn/dMi over M^, < Mi < 00 
to obtain I^pbh = e^nh x^'^'dx = 2 x 10-^£^{nh/0.l), where x = Mi/M^. Finally, 
assuming Qh = 0.1 Eq. (|D leads to an U.L. of QpBu < 6 x 10~^, being comparable to 
that from the diffuse 7-ray background spectrum 10^^) 0. 

In closing, if future long-duration (~ 8 days) balloon flights allow us to precisely measure 
the cosmic-ray p flux at kinetic energies from 0.2 to 2 GeV, the U.L. on TZ can be significantly 
reduced to ~ 2 x 10^^ pc^'^yr^^, or more importantly, such observations could confirm the 
existence of evaporating PBHs. 
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FIGURES 

FIG. 1. (a) Distribution of Dq{z)~^ in halo models I-III as shown in the z-direction where the 
Solar system is located on the galactic plane (z = 0). (b) Distribution of PBHs contributing to 
the local p flux F near the Solar system, obtained by assuming = 1.0 x 10^*^, as shown in the 
z-direction. (c) The same as (b) shown in the r-direction where the Solar system is located at 
r = 8.5 kpc. 

FIG. 2. (a) Simulated local interstellar PBH-p flux using model III with TZ = 1 x 10^^ and 
1 X 10^^ pc^^yr^^. (b) Corresponding p/p ratio for PBH-p's at the top of the atmosphere (dashed 
lines), the expected p/p ratio for secondary p's (solid line) IQ, and observational data p|,p^|-p3|. 
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